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ABSTRACT 

Recent astrophysical observations have provided strong evidence that the 
present expansion of the universe is accelerating, powered by the energy den- 
sity associated with a cosmological term. Assuming the latter to be not simply 
a constant term but a "quintessence" field, we study the radiation of quanta 
of such a "quintessence" field ("quintons") by binary systems of different types 
and compare intensities to those of standard tensor gravitational wave emission. 
We consider both the case in which the quintessence field varies only over cos- 
mological distances and the case in which it is modified spatially by (strong) 
gravitational fields, a condition that results in bounds on the gradient of the 
scalar field. We show that, in both the first case and, because of a bound we 
derive from the Hulse- Taylor pulsar, in the second, there is not sufficient quinton 
radiation to affect expected LISA and LIGO gravity wave signals from binary 
systems. We show that, in the second case, the Large Hadron Collider is capable 
of setting a bound similar to that from the binary pulsar. 

Subject headings: radiation mechanisms: general — gravitational waves — 
binary systems 



1. Introduction 

The past few years have witnessed a revolution in cosmology, a field that has been 
growing from data poor to data rich by leaps and bounds since the early 80's. This revolution 
followed the discovery, by the study of distant supernovae (Garnavich et al. 1998, Perlmutter 
et al. 1998), that the expansion of the universe is currently accelerating. A similar conclusion 
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was reached in the "concordance analysis" of the cosmic microwave background (CMB) data 
by the WMAP collaboration (Bennett et al. 2003; Page et al. 2003) which concluded that 
the universe has a value of fl, the ratio of its density to the critical density, consistent with 
Q = 1. In addition, the analysis of the WMAP data allowed an independent estimate of 
the contribution to this value by matter alone (including that of the "dark matter" which 
makes up most of the gravitationally bound matter); this was found not to exceed 30%, 
suggesting that a cosmological constant term could be responsible for the remainder, leading 
to an estimate Qa ~ 0.7. 

The presence of a cosmological term of this magnitude has been the cause for great 
concern in the field: The "natural" value of such a term in the context of theories of the 
fundamental interactions would be huge, of order of the Planck energy density, the only 
"typical" value in a theory that contains a mass scale such as the Planck mass Mp (Carroll 
2000; Peebles & Ratra 2003). The hope, then, was that a (hitherto unknown) conservation 
law would set its value to precisely zero, a value generally consistent with the pre-1998 
data. The presence of a small, non-zero value for the cosmological term is then seen as 
another "fine tunning" problem in cosmology. Nevertheless, independently of the issue of 
the magnitude of the cosmological term, it is generally expected that this term does not 
represent a universal constant associated with the Lagrangian of gravitational interaction. 
Rather, it is thought that it hkely varies like a field, whose value is constant on cosmological 
scales, while its magnitude is varying slowly in time. This field, being quite distinct from 
the other known fields and forms of matter has been given the name quintessence (Caldwell, 
Dave & Steinhardt 1998), borrowing the nomenclature from Aristotle's "fifth substance" 
that was supposedly involved in the make up of the universe. 

Once, however, one decides that the term that drives the acceleration of the universe is 
not an a priori constant but a (scalar) field, one is immediately led to the notion that such 
a field obeys its own (scalar wave) equation and that it should be neither constant in time 
(a feature generally used in the literature) nor uniform on non-cosmological scales. As such 
it carries its own "kinetic" energy and potential energy which couple to the gravitational 
field thereby affecting both the metric around a bound object and also being affected by it, 
through the contributions of the metric to the covariant derivative. 

Interestingly enough, an exact solution of the combined set of Einstein and zero-mass 
scalar field equations in the static spherically symmetric case has been derived independently 
several times in the literature, (e.g. Buchdahl 1959) including by one of us (Mannheim & 
Kazanas 1991 ). This last reference in particular, conjectured that the scalar field considered 
is none other than the Higgs field of high energy physics, which presumably is responsible for 
giving masses to particles through a spontaneous symmetry breaking process. This is rather 
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relevant in that our present study involves also the study of its effects in high energy collisions 
(of course, this field is expected to make the usual (huge) contribution to the cosmological 
constant, which must be somehow cancelled as discussed above; the effects discussed by 
Mannheim & Kazanas (1991) pertain only to the effects of the kinetic energy of the scalar 
field). 

Motivated by the above considerations, we examine, in the present note, the case in 
which both the scalar and gravitational fields are space and time dependent, i.e. the case 
of quinton emission by their combined action. Such a study becomes imperative with the 
development of detectors of gravitational radiation that either currently arc (TAMA) or 
about to become operational (LIGO), or are planned to be built in the not too distant future 
(LISA). Identifying gravitational waves in these facilities will depend on the comparison of 
the data to theoretically-derived, computer-generated, wave forms. These forms, and hence 
the interpretation of the data, could be impacted by competing processes not accounted for 
in template generation. 

Due to the nature of our study, we have proceeded classically in the linearized regime 
with our main goal seeking conditions and physical situations for which radiation of scalars 
might become comparable to that of gravitational waves. To this end we explore both 

astrophysical systems, as well as systems in the laboratory (e.g. high energy collisions). 
We will also explore both the case in which the variation of the scalar field is limited to 
cosmological times and distances and the case in which the scalar field configuration is 
infiuenced by the presence of nearby (compact) objects. 

Since we will be comparing scalar radiation with graviton radiation, we recall here 
(Shapiro & Teukolsky 1983) classical order of magnitude formulae for the latter. For a binary 
system, gravitational radiation luminosity is given in terms of the reduced and total masses, 
Ij, and M respectively, by the expression 

Law = (32/5)G^MVV(c'a^) (1) 

with M — Ml + M2, 11"^ — Mf -I- M^-*^, and a — ai + a2 where Mj and are the mass 
and distance from the center of mass (and origin) of body i {i — 1,2) The rate at which the 
period, u'^ — P — 27r/a; {co^a^ — GM), changes is given by 

p-^dP/dt = -(96G'V5c^)MV/a^ (2) 

In §2 we set out the relevant equations and their linearized form while in §3 we ex- 
amine the intensity of scalar field radiation resulting from the changing gravitational field 
(or metric) of a binary system and compare it to that of the gravitational radiation. We 
work, first, within the approximation that the (cosmological) scalar field is not affected by 
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the gravitational field of the binary system. In §4 we examine the competition between 
gravitational and scalar radiation taking into account the change in the scalar field due to 
the gravitational field of the binary using the solution of Mannheim & Kazanas (1991) and 
finally, in §5 we present and discuss our conclusions. 



where the box operator denotes the covariant d' Alembertian, F'^ are the Christoffel symbols 
associated with the (in general) curved space into which ^ operates and m is the (effective) 
mass of the (quinton) field ^. Our calculations are based on the approximation of first 
solving Equation (3) for the "semi-static" case and then computing the radiation field by 
approximating the right hand side by zero and using the F-term, with the semi-static quinton 
field, as the quinton radiation source. 

Assuming space-time to be sufficiently close to flat, one can use a perturbative approach 
to calculate the components of the metric tenson g^^'^ and then the Christoffel symbols. Thus, 
in the slow velocity regime, one can expand the metric tensor to powers of v/c, where v is 
the magnitude of the velocity of the matter components involved. In this approximation, to 
second order in f/c in g^^ and to first order in f/c in (7**, the space part of the metric tensor, 
the departure from fiat is equal to twice the gravitational potential (p. In this approximation 
the diagonal metric tensor components and the square root of its determinant are (Weinberg 
1972), remebering that —g — det{g^i,), 



5l = [(-l + 20),(l + 20)], and (-5)^/^ = [(1 - 20)^] = 1 - 20 + 0(0^) (4) 



2. Basic Equations 



Let denote a scalar field obeying the scalar field equation 



(3) 



leading to 



[-l + 40 + O(02),i + O(02)] 



(5) 



which then yields for the Christoffel symbols 



(6) 



With the above expression for F'*, the equation obeyed by the scalar field in curved space- 
time becomes 

{d'^ - V^)* - 4 9o05o^ = -m^* ~ (7) 
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This equation is a wave equation with an inhomogeneous term. Its solution can be 
found in standard texts (e.g. Jackson 1962), i.e. 

t) = 4 [ -J^^doM^', t')do<P{x', t') (8) 
\ X — x' \ 

where is the unperturbed (cosmologicany varying only) scalar field at point x' and time 
t' and we are suppressing factors of c except when evaluating expressions numerically. 

Considering the combined gravitational - scalar field variations appropriate for a binary 
system in circular orbit centered on the center of mass of the two bodies, the above equation 
reduces to 



X 



\x' — x\ 



— * — * 

\x' - ai[ti)\-^ F'-a2(^2)r 



(9) 



t' ^t-\x-x'\^t-x + x-x' (10) 



*«(x,i)= 4^(50*0) 
where 
and 

ti = t' - \x' - ai{ti) \ c^t' -x' + x' ■ ttiiti) (11) 

where (liiti) is the "doubly retarded" vector connecting the center of mass of the system with 
body i. In other words ti is the time at which a light signal from aj(ij) would have to leave in 
order to get to x at time t via x' at time t —t — \x — x'\. We will also use do^ — mefr^o, with 
mlo^l ^ Ude: where Ude — lO^eV /cm^ is the observed dark energy density The quinton 
(scalar field) at x is the result of coherent addition of the ^ generated by the time dependent 
field at x' from the masses at ai{ti). 

In the long wavelength approximation {u!~^ > a), most of the contribution to the integral 
comes from distances less than x' ~ Tr/cu.We take |x — x'| ~ x for the outer denominator. 
We use 

ai iU) =ai {t' - \x'- ai {U)\) ~ a, f.i'-i.t' -a,{u)\) ^ - + ^^^/. (12) 

In the same (long wavelength) approximation, we can ignore differences between the time 
dependence of ai and 02 beyond that included in Equation (12). We now perform the time 
derivatives of the a' s in Equation (9) making use of Equation (12). The result is 

■^r{x, t) = 4G'C/J/Je'^'"(*-^)x-' j dx'dn'iu'' i' [{i'- aifMi + {i'- ^2)'M2]} (13) 



Other terms, at least in lowest order approximation, vanish and/or have the wrong time 
dependence. We make the approximations: (1) that the lower limits in the two x' integrals. 
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ai and 02 are both approximately a/2, with a — ai + a2 and (2) that the upper hmits are 
of order n/uj past which the oscillating exponentials dampen any contribution. We use the 
fact that Miof + M2a| = /la^ to obtain 

^r{x, t) = (167rV3a;)e2^'*^(*-^) G U]ila''uJiJi (14) 

With this result, the fact that the intensity L^, of quinton radiation is Attx'^T^^, and the 
expression for the Oi-component of the scalar field stress-energy tensor — df ^ rBj.'^ r, we 
obtain for the intensity 

Lvp ~ lO^f/DsCVaV (15) 

where we have integrated over 9 and averaged over This expression is applied to several 
specific cases in the next section. And, in the section after that, the analogue to this result, 
for the case in which the scalar field solution is modified by the gravitational field of the 
compact system, is derived and applied. 



3. Slowly Varying Quintessence Field 

We consider the case in which the scalar (quintessence) field (in section 11) is not modified 
by the presence of strong gravitational fields. In the following section we consider the case 
studied by Mannheim & Kazanas (1991) in which the scalar field is modified by changes in 
the gravitational field. 

We begin by applying Eq. (15) to a binary system in circular motion of radius a and 
Keplerian angular frequency uo. Using Kepler's law uj"^ = GM we obtain in terms of 
Ms = M/MQ,fis = I^/Mq, and aw = a/10^° cm 

= 8.7 X 10^UDEG'^M^iJ,^/(c'^a'^) ~ 10^(//, M,/aio)^ erg/s (16) 

where we have inserted needed factors of c in Eq. (15) [{Uj:)ec){G'^/c^)]. Comparing with 
Eq.(l), Lg ~ 1.7 X 10^^ M^ii'l/alo erg/s, we see that 

Thus quinton emission will dominate emission of gravitational radiation only for very large 
values of the orbit radius at which point they are both negligible. 

Looking at larger mass objects, i.e. galaxies, clusters, superclusters etc... we can set 
M ~ (4/3)7ra3pcQM or ~ lO-^^a^Q. Eq. (17) requires > lO^^M^^ for quinton 
production to win. For two galaxies of M ~ this gives 10^ hght years, somewhat 
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larger that the average galaxy separation. For a star in a galaxy circling the center of mass, 
is negligible compared to Lq- 

Going the other way, we can look at small systems, asking how much energy is radiated 
in quintons in an excited state lifetime of an atom. Going back to Eq. (15) we insert 
uja — V — ac/uB, where ub is the principal quantum number. Taking /i to be the mass of 
the electron, we see that 

= 10"^^^ eV/s (18) 

We can also ask for the enchancement that would result in the case that there are large 
compact dimensions (Arkani-Hamed, Dimopoulos & Dvali et al., 1998) and the quintons can 
travel in the bulk. Following that reference for the case n = 2 of just two extra compact 
dimensions, we would get from Eq. (9), for an atom (a ~ 10~^), 

The result would be to modify Eq.(18) by a factor of ~ 10^^ which does not appear enough 
to make it of experimental interest. There is a larger effect if we consider a smaller (nuclear) 
system such as a— decay of a long-lived isotope. In that case, a 10 MeV a— particle could 
give (a;a)^ ~ 10~^c^ so that Eq. (16) would be 

Ls, = 10\lO'eV/ cm')c ( ( -^V ("^^V (lmm/5 • lO-i^^m)^ ~ 10"^^ eV/s (20) 

\mpij \mpij \ c J 

Thus, even a mole (10^^) of an isotope with a billion year half life would have less than an 
electron volt of energy loss into quinton radiation. Finally, we turn to accelerator production, 
p — p collisions at the LHC. We take a ~ cr^^^ ~ 10~^^ cm and compute the energy radiated 
into quintons in a collision. 

A£; = L^5t (21) 

with 5t — a/c and 



= 8.7 X lQ^{UDEG'^){ix'^/h){aw/c)\lmm/aY (22) 

where we have, again, assumed two extra compact dimensions. We have also assumed 
auj = c, although it is possible that, in a quantum treatment, we might have au — > aE 
giving a much larger result. We have, of course, set ji = E. Again the result is small: 
/S.E ~ lO"^^^^^. Thus, based on a quinton field varying only over cosmological times, 
quinton radiation does not approach gravitational radiation for any of the 3 cases considered 
in atomic and nuclear transitions and high energy collisions, nor in the astrophysical binary 
systems considered. 
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4. Quintessence Field Varying in Strong Gravitational Field 

We turn now to the possibility that the scalar field is modified in the presence of grav- 
itational fields. Specifically, we address the cases of Section III in light of the results of the 
work of Mannheim & Kazanas (1991). These authors considered Eq. (3) written in the form 

(3^9„[(-s)"^a"*l = %^0 (23) 
along with Einstein's equations written in the form 

1 _ Qt^^Q<^^ + -g^^'d^^Ifda'^ = -g^^Vi"^) (24) 



SttG T 

They find, as did Buchdahl (1959), closed form solutions for the case V(^') = 0. With 
V{^) ~ r7igg\['^ and m^<s ~ 10"^^ eV, as demanded by the dark energy observations, their 
two solutions should be good approximations. One is that ^ is a constant, unaffected by 
gravitational fields. This, of course, is just the case considered above. The second solution, 
of the coupled equations, is 

ds^ = -H{p)dt^ + J(p) [cip2 + p'^dn] (25) 
The functions in the above equation are 

*(p) = ^/n ( ^ ~ ^° ) + constant (26) 



2ro VP + ^o 
= (^) 

2 / X -d/2ro 



where 

d = 2MG = 4{rl - TiGK'^f'^ (29) 
and it can be shown (Mannheim & Kazanas 1991) that Tq is restricted to the region 

MG/2 <ro< MG, (30) 

indicating that tq is of the same order of magnitude as the Schwarzschild radius. We set 
ttK'^ = GM^(3^, ro = MG-f with < /3 < 1 and 1/2 < 7 < 1. We assume that the 
dimensionless quantities, f3 and 7, are independent of the gravitational field, that is, like 
Newton's constant, are the same for all masses. 
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Using the above equations we compute F^, ^^h, L$ and L^/Lqw- We first obtain 

3-"%(/'V) = r- (31) 

= 4fl-4)7^)"'''(P^-6'-?--'-oP) (32) 
2 

~ - for p > a » To (33) 
P 

Similarly, we have = (2/p)9op, ~ K/p^, and ~ {K/p^)dop. 
We rewrite the solution of the wave equation (3) as 

•^mk{x, t) = x-^ j dVr'^(5, i')9^*o (34) 

for each of the two bodies in the binary system separately. Using P—x' — ai{ti), we will add 
the two contributions to "^mk- The time and p components give much different results for 
"^MK and Lmk- For the time component, do P is simply do a and the leading contribution 
is 

*MK,o = x-^ J d^x'i'iiuj ai ■x'){2iuj ■x')Kjx'^ = {"i^n /?,x)e^''^^*-''^ K^aluHn{l/ aw)+{l 2) 

(35) 

which gives, after using T:{Kia\ + K^a^Y — l3'^Gp?a^, 

Lmk,q = {6i7r/3)^p\Gp^){auj)*ujHn\auj/c) (36) 
The p contribution is given by 

*MK,p = {2K,/x) J dV/l x' - ^1 |3 + (1 ^ 2) (37) 

We expand the denominator(s) in powers of the a^; noting that the term of first order in a, 
vanishes after the angular integration, we are left with terms of order 2. The result is 

'^MK,p = (6Xi/x) J d^x'[A{i'- ai)^ - al]/x" + (1^2) = 27r(Ki + K^) (38) 

from which we get 

Lmk,p = 2(47r)2G/x2a;2/?2 (39) 
Here, we have approximated 7r(i^^ + Xf) by its equal mass value. 

Clearly, the relative factor of (acu/c)^ between Lmk,o and Lmk,p means that the latter 
will be more important for astrophysical objects and as well as atomic and nuclear ones. 



while the former will dominate in high energy physics. Turning first to astrophysics, one 
may note 




5/?V4 



(40) 



where v is the velocity of the lighter member in the non-equal mass case. Wc note the 
similarity between Eq.(39) and Lqw- This might be expected since there is no additional 
mass scale, such as the t/o^; in section II which enters in forming "i/ji. . 

We can now estimate an upper bound on /3 from the Hulse - Taylor pulsar which has 
an orbital period of 28,000 sec and the fact that the rate change of its orbital period agrees 
to better than one percent with the prediction of General Relativity. From the fact that 
Lmk/Lgw < 0.01 wc obtain that < 2 x 10"^(f/c)*^ ~ 10~^^. While this appears a 
stringent bound on the constant K, understanding its full implications awaits simultaneous 
solution of Equations (23, 24) in the interior region as well as the exterior region which 
should permit evaluating K and Tq in terms of the interior mass distribution. 

Turning to the high energy case, wc compute the energy that would be lost into quinton 
radiation in proton-proton collisions. Lmk,o of Equation (36) above dominates. In it, we 
take // to be 7 TeV as with the LHC in section 3 and cu to be c/a, while recognizing that 
it might be larger in a quantum treatment. Again, we multiply L by a/c to obtain the 
energy radiated during the encounter. We include enhancement from two extra compact 
dimensions. We ignore the logarithm. These give 



This imphes that, if significant unexplained energy loss in p — p collisions at the LHC is not 

found, a bound on (3"^ slightly better than that above from the Hulse- Taylor pulsar can be 
inferred for models with compact extra dimensions. It should also be possible to infer such 
a bound from cosmic ray data from the Auger project (www.auger.org) which will study 
p-p scattering at about 30 times LHC energies (at the center of Mmss). One might even 
infer from the existence of ultra high energy cosmic ray observations that proton-proton 
interactions do not lose large amounts of energy into unobserved particles, implying a bound 
on of about 10~^^. On the other hand, the cosmic ray spectrum does exhibit a feature 
known as the "knee" at an energy ~ 10^^-^ eV which corresponds to roughly 1 TeV 
at the center of mass. This is a steepening in the slope cosmic ray spectrum by 0.3 - 0.4 
over half a decade in energy. The very limited energy range over which this change in the 
cosmic ray spectral index occurs essentially precludes its explanation as simply a cosmic ray 
propagation effect. In fact Kazanas & Nicolaidis (2003) suggested that this feature heralds 



AE = {6ATi/3f{E/mpifi3\uj^a/c){lmm/ay = lO'^'P^GeV 



(41) 
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the emergence of physics beyond the Standard Model. In particular, these authors have 
argued that such a feature is the result of energy lost in cosmic ray collisions to gravitons, as 
suggested by the theories of extra large dimensions which presumably have the graviational 
constant increase with energy to that of strong interactions at energies ~ 1 TeV. It would not 
be in conflict with our Hulse- Taylor pulsar bound, if quinton emission were also to contribute 
to this cosmic ray feature. 

5. Summary 

General relativity requires that any quintessence field couple, through the covariant 
derivative, to the gravitational field. Thus, any time varying gravitational field will produce 
quintons. The rate of production, however, will be proportional to the square of the derivative 
of the quintessence field. We have evaluated that rate, in a classical approximation, in Section 
II, for the case in which the field only varies over cosmological times. We apphed the results 
to a variety of binary systems in Section III, but found no cases of interest in which energy 
loss through quintons dominated energy loss through gravitons. 

In Section IV we turned to the perhaps more realistic case in which the space variation 
of the quintessence field is affected by the presence of a gravitational field. We worked 
with the exact (exterior) solution for the massless case as written down by Mannheim and 
Kazanas (1991). The result, in that case, was quinton production dominated by different 
components for low velocities of the binary system than for high. We were able to bound 
an integration constant in the solution cited by requiring that quinton emission from the 
Hulse- Taylor binary pulsar system represents less than one percent of the energy loss. That 
limit on the parameter was sufficient to make it difficult to identify observable effects in the 
astrophysical phenomena considered. We were also able to show that, with the assumption 
of large, compact extra dimensions similar and stronger bounds could be derived from data, 
when available, from the LHC and project Auger (and perhaps from current cosmic ray data 
showing the existence of ultra high energy cosmic rays). 

In summary, it appears that quinton emission is unlikely to be of importance in inter- 
preting signals received by LIGO or LISA, or in laboratory experiments if the quinton field 
has no coupling to matter beyond the indirect coupling provided by the covariant derivative 
or if large compact extra dimensions do not exist. 

It is a pleasure to acknowledge helpful conversations with Richard Fahey, Breno Im- 
biriba, Doris Rosenbaum and Michael Turner. 
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